Abrikosov vortex escape from a columnar defect 
as a topological electronic transition in vortex core 
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We study microscopic scenario of vortex escape from a columnar defect under the influence of a 
transport current. For defect radii smaller than the superconducting coherence length the depinning 
process is shown to be a consequence of two subsequent topological electronic transitions in a trapped 
vortex core. The first transition at a critical current jl is associated with the opening of Fermi surface 
segments corresponding to the creation of a vortex-antivortex pair bound to the defect. The second 
transition at a certain current jd > Jl is caused by merging of different Fermi surface segments, 
which accompanies the formation of a freely moving vortex. 

PACS numbers: 



The study of magnetic and transport properties of 
type-II superconductors in the presence of artificial pin- 
ning centers is known to be an important direction in 
the physics of vortex matter [1| . An obvious first step to 
understanding of the pinning related phenomena is asso- 
ciated with the evaluation of the individual pinning force 
acting on a vortex in the presence of a defect. Various 
approaches used for such analysis are based on the obser- 
vation that a defect attracts the vortex in order to avoid 
loss of condensation energy in the core and to decrease a 
kinetic energy due to perturbations of screening currents 
and magnetic field induced by the defect 0-01 ■ Close to 
the critical temperature Tc both these effects can be de- 
scribed within the Ginzburg-Landau (GL) theory valid 
for large length scales well exceeding the superconducting 
coherence length ^ at zero temperature. 

Such long wavelength approach is, of course, no more 
valid for low temperatures T Tc and small defect 
size a < ^. Microscopic approach for calculation of the 
pinning energy and appropriate modification of the GL 
functional for a particular case of a small point impu- 
rity with the scattering cross section CTsc ^ S,^ have been 
previously studied in Ref . 0, 0] . Considering the mi- 
croscopic theory one should take into account the be- 
havior of the subgap fermionic states bound to the vor- 
tex core which are known to determine both the struc- 
ture and dynamics of vortex lines in the low tempera- 
ture limit. These subgap states are known to form the 
so-called anomalous spectral branch crossing the Fermi 
level. For well separated vortices the behavior of the 
anomalous branches can be described by the Garoli-de 
Gennes-Matricon (CdGM) theory Q: for each individ- 
ual vortex the energy eo(M) of subgap states varies from 
— Aq to -f Aq as one changes the angular momentum fi 
defined with respect to the vortex axis. Here Aq is the 
superconducting gap value far from the vortex axis. At 
small energies \e\ ^ Aq the spectrum is a linear func- 
tion of fi: eo{fi,k±) ~ —fxhoj, where Huj « AolnA/(fcj_^) 
is the interlevel spacing, ^ = hVp/'^o, k\ kp — k^, 
kz is the momentum projection on the vortex axis, kp 



and Vf are the Fermi momentum and velocity, respec- 
tively. The logarithmic factor with A ^ Aq/T appears 
at low temperatures due to the so-called Kramer-Pesch 
effect [l^ . Neglecting for T ^ fiw the quantization of 
the angular momentum fi one can consider the anoma- 
lous branch to cross the Fermi level at = for all 
orientations of the Fermi momentum hp. Thus, in the 
space iJ, — \ip we obtain a Fermi surface (FS) for excita- 
tions localized within the vortex core (see for review). 
Ghanging magnetic field or transport current we can get 
switching between different vortex states which is accom- 
panied by the changes in the FS topology. Such topolog- 
ical transitions in quasiparticle spectra of vortex systems 
are similar to the well-known Lifshits transitions which 
occur in the band spectra of metals jj, |l3[. One can 



separate two generic examples of such transitions in vor- 
tex matter: (i) opening of FS segments corresponding to 
the creation of vortex-antivortex pairs 14, 1^; (ii) merg- 
ing and reconnection of different FS segments via the 
Landau-Zener tunneling 15,3]. The basic properties of 
vortex matter such as pinning and transport characteris- 
tics, heat transport in the vortex state and peculiarities of 
the local density of states should be strongly affected by 
these changes in the FS topology. The goal of this Letter 
is to suggest a theory of topological electronic transitions 
which occur in a pinned vortex core under the influence 
of transport current. 

The mechanism of these transitions is closely related to 
the effect of elastic scattering on the CdGM levels. Re- 
sulting modification of the anomalous spectral branch is 
noticeable even for impurity atoms introduced in a vortex 
core [TtI and becomes much more pronounced provided 
we consider defects of the size well exceeding the Fermi 
wavelength. The absence of the FS is a hallmark of a 
vortex pinned by an insulating columnar defect and it is 
natural to expect that vortex depinning should be accom- 
panied by the FS formation. For a vortex line trapped by 
such defect the subgap spectrum has been analyzed re- 
cently within the quasiclassical approach in Ref. [18| . For 
rather large angular momenta > iia = k±a quasiclas- 



2 





'J 


V . 






a/ 




a) 









A 




m 


8 

a 




-A 

)---°- — 


E 

a 


■ 1 ■ 1 ■ 

A 

m 







Figure 1: Fig.l Different types of quasiclassical trajectories 
(a) and spectrum of subgap states (b) for a vortex pinned by 
tlie defect. 

sical trajectories of electron-hole excitations do not expe- 
rience reflection at the defect surface and, thus, the sub- 
gap spectrum coincides with the CdGM one. In the op- 
posite case 1^1 < [la the normal reflection of trajectories 
at the defect destroys the low energy part of CdGM spec- 
tral branch: with the decrease in below the threshold 
value \ia the energy of the subgap bound state rapidly ap- 
proaches the superconducting gap Aq (see Fig. 1). Thus, 
contrary to the CdGM case the subgap spectral branch 
of the pinned vortex does not cross the Fermi level: there 
appears a minigap Am(a) ^ |eo(fcFa, ^f)! in the quasi- 
particle spectrum. The spectral branch Ea for trajectories 
touching the defect can be approximated by vertical lines 
passing from ±Am to ±Ao for [i = =F/^q (see Fig. lb). 

In this Letter we propose a microscopic scenario of 
vortex escape from a columnar defect of a radius a smaller 
than the coherence length ^ under the influence of an 
external supercurrent j = engV applied perpendicular 
to the vortex axis (Fig.l). The transport supercurrent 
results in a Doppler shift = (fi.k-V) of the quasiparticle 
energy and the subgap spectrum takes the form: 

e(/i, k±,9p) ~ e{fi, k±) + hk±V cosdp , (1) 

where e = — Aosign/i for quasiparticles with \fi\ < fia 
and e — So for quasiparticles with > jia- Here the 
quasiparticle momentum k = k±{cos6p, sin^p) defines 
the trajectory orientation angle 9p in the {x, y) plane 
and V = yxo. The depinning process can be viewed as 



two subsequent topological electronic transitions in the 
vortex core which restore the anomalous spectral branch 
crossing the Fermi level and repair the FS. The first tran- 
sition at the critical current jl = sUsVl is associated 
with the opening of separate segments of the FS which 
appear when the Doppler-shifted branch of the quasipar- 
ticle spectrum crosses the Fermi level and the minigap in 
the quasiparticle spectrum vanishes: 

V>Vl = Amia)/hkF . (2) 

The latter condition coincides, in fact, with the famous 
Landau criterion of superfiuidity [l^ . The corresponding 
critical current density can be expressed via the depairing 
current density jc — ekpAo/^n^h and CdGM spectrum: 

Jl = jc leoikpa, fcF)|/Ao < jc ■ (3) 

Considering the zero temperature limit and assuming 
a <^ we obtain the following expressions for the critical 
velocity and current: 

Vl - Aoa\nA/hkF^ , - Jc(a/0 In A < , (4) 

where A ~ ^/a. Taking e = in Eq.(IT]) wc find the 
equation defining FS segments in the (p, Op) plane: 

£Q{n,k±_) = -hk±y coaOp , > /ia . (5) 

These segments are shown in Fig. 2 by solid lines in the 
limit a ^ ^ when £o(Mi ^-l) = —huj^. The ends of the seg- 
ments are joined by the lines ^ — ±fia which correspond 
to the spectral branch e^. The direction of quasiparticle 
trajectory precession along the resulting closed classical 
orbits in the jj, — 9p plane is determined by the Hamil- 
ton equation: hdOp/dt = de/dfji. The formation of FS 
segments is obviously associated with the appearance of 
superconducting phase singularities outside the columnar 
defect. 

We now proceed with the analysis of the structure of 
a resulting vortex state bound to the defect. For the 
sake of simplicity hereafter we put k± — kp which is ap- 
propriate for strongly anisotropic layered material with 
vortex lines oriented along the anisotropy axis. One can 
separate two contributions to the gap function A(r, 9) in 
the self-consistency equation: (i) the term Ai associated 
with the subgap states localized within the vortex core; 
(ii) the term A2 which is defined by the quasiparticle 
states with large energies e > Aq. We introduce here the 
polar coordinate system (r, 9) in the plane perpendicular 
to the vortex axis and the origin chosen at the cylindrical 
defect center (see Fig. 1). At rather low temperatures the 
term A2 is only weakly affected by the transport current 
j ^ jc since the corresponding Doppler shift is much 
smaller than the superconducting gap Aq. On the con- 
trary, the first contribution to the gap is strongly modi- 
fied in the presence of a transport current provided the 
Doppler shift exceeds the minigap Am(a)- Our further 
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Figure 2: Fig. 2 Schematic plot of Fermi smface segments 
in the {fi,9-p) plane for Vl < V <^ Vd (solid blue line) and 
V > Vd (dashed red line). 



calculations of the current induced correction to the term 
Ai{r,d) are based on a semiclassical version of the self- 
consistency equation (see [2C 



<5Ai(r, 



27r 



e'^" (sign e — sign e) , (6) 



where £(/i) {e{fi,0p)) is the subgap spectrum in the ab- 
sence (presence) of transport current, /i — kpr sm{6p~6), 
and Af, ^ Aq is the amplitude of the gap value contribu- 
tion coming from the CdGM spectral branch. Here we ne- 
glect the coordinate dependence of normalized electron- 
and hole-like wavefunctions at the scale r ~ a ^ ^ and 
consider the zero temperature limit. Close to the thresh- 
old velocity Vl the spectrum e can be considered pertur- 
batively in the form ([T]) with k± = kp. Note that here we 
restrict ourselves to the small /i limit assuming j <C jc- 
For small superfluid velocities V <Vl, the gap function 
Ai = Af,e'* is undisturbed {SAi = 0). Otherwise, if 
V > Vl, the Doppler shift results in qualitative changes 
in the gap function distribution. A resulting piecewise 
regular gap function 



Ai(r, 



|Ai|exp(i0) = Afce''' + ,5Ai(r,0) 



is defined in a set of regions shown in Fig. 3: 



Ai(r,0) = Afac'' 



where sin ar 



1, 

1 - 2Vr2 - a^/r, 
1 - e*"- + e-*"o , 



in A or C , 
mB, 
in D, 
mE, 



(7) 



sm ao 



rg cos 6 



y^r'^ + 2ror sin + Tq 



and ro = V/u. By dashed lines in Fig. 3 we also 
show schematically the distribution of the superconduct- 
ing phase gradient V0 around the defect. One can clearly 




Figure 3: Fig. 3 Schematic distribution of the superconduct- 
ing phase gradient (dashed lines) for a bound vortex con- 
figuration. 



observe strong changes in topology of the phase (f> distri- 
bution caused by the superflow: there appear a phase 
singularity at the point S {x = 0, y = —tq) and the line 
L {r = 2a /V3) of zeros of the superconducting order pa- 
rameter Ai in the region B (see Fig. 3). The singularity 
at the point S has the same vorticity as the initial vor- 
tex captured by the defect. With an increase in the su- 
perfluid velocity V this singularity moves away from the 
defect. One can see that the modification of the phase cj) 
distribution shown in Fig. 3 corresponds to the formation 
of a bound vortex- antivortex pair. Note that the vortex 
and antivortex positions found here coincide with the po- 
sitions of the phase singularities of a total gap function 
since the amplitude of the second part of the gap func- 
tion A2 is rather small: |A2(r ^ ro)| ^ AqTo/^ ^ A;,. 

Certainly well above the threshold velocity Vl the gap 
profile and stable vortex configuration should be deter- 
mined from the self-consistent procedure taking account 
of contribution of delocalized excitations. It is natural to 
expect that the instability of the bound vortex state and 
appearance of a freely moving vortex should be accompa- 
nied by the formation of a continuous anomalous branch 
connecting the energies below — Aq and above -l-Ag. It is 
this branch which provides necessary conditions for spec- 
tral flow through the energy gap and, thus, for energy dis- 
sipation during the vortex escape from the defect. The 
obstacles on the way of quasiparticle flow through the gap 
could be removed provided the FS segments would merge 
and reconnect to form a cosine curve fj-{0p) typical for a 
free vortex with a continuous anomalous spectral branch 
(compare solid and dash lines in Fig. 2). This transforma- 
tion causes topological changes in FS geometry because 
new cosine FS is open and can not be contracted to a 
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point by continuous transformations. Correspondingly, 
one can suggest three obvious scenarios of destruction of 
the bound vortex configuration: (i) first, the FS segments 
can merge forming an open FS due to quantum mechan- 
ical tunneling; (ii) second, the spectral fiow through the 
energy gap can occur due to the impurity scattering be- 
tween different Fermi surface segments; (iii) third, the 
free vortex can be formed when the Doppler shift reaches 
the gap value, i.e., the current density approaches the de- 
pairing one. The first scenario can be realized due to the 
Landau-Zener tunneling between classical orbits and can 
be understood within the quantum mechanical picture 
based on the commutation relation for canonically con- 
jugated variables n and 9p: [fj.,9p\ = —i. Here we use 
square brackets to denote a commutator of two opera- 
tors. The efficiency of this tunneling is determined by 
the ratio of the distance 2/iQ between the orbits to the 
angular momentum uncertainty A/x. The latter value can 
be estimated from the uncertainty principle AfiAOp ^ 1 
and the equation = k±V cos9p/uj describing the classi- 
cal orbit for small /x values. Considering the trajectories 
with 9p close to ±7r/2 we find: A/i = k±VA9p/Ld and 
A/i ^k\Y jijj. Thus, Landau-Zener tunneling results 
in the merging of two FS segments provided the condition 
Ma < \/k±y/(jj is fulfilled. This second transition in the 
spectrum accompanied by the change in the Fermi sur- 
face topology occurs at the critical velocity and current 
given by the relations: 

Vd = kpuja^ , id = jcikpo^/O In A = kpajL > Jl ■ (8) 

For j > id the FS is described by the Eq. [5] for all 
/X values (see dashed line in Fig. 2) and corresponds 
to a Doppler-shifted spectrum of a free vortex: e — 
—f)w^ + hkpVco&9p for |/x| ^ kpS,. The continuous path 
for the quasiparticle spectral fiow from the energies be- 
low — Aq to the ones above -I-Aq is now restored providing 
conditions for a dissipative vortex motion. 

Second scenario of the recovering such continuous path 
through the gap is associated with impurity scattering. 
The scattering rate I/ts between the quasiparticle states 
at different FS segments shown in Fig. 2 should be re- 
duced compared to the impurity scattering rate l/r in 
a normal metal due to the small size 69p of FS seg- 
ments in the phase plane fi — 9p: 1/ts ~ S9p/T, where 
cos{S9p/2) ~ uja/V. Provided the scattering rate 1/ts 
becomes comparable with the minigap A^ = kpahu 
for a pinned vortex the FS segments can be no more 
considered as isolated ones. This condition imposes 
a restriction on the size of the isolated FS segments 
156*^ < 59max = lAmT/h < TT and gives us an estimate 
for the upper critical velocity and current which destroy 
the bound vortex state: 



COs{59max/'2') ' COs((50„ax /2) 

where the constant 7 is of the order unity. Note that 




Figure 4; Fig. 4 Schematic dependencies of critical current 
densities vs the defect radius a for clean (a) and dirty (b) 
limits. 

in the limit ^^A^Tjh > tt the impurity scattering effect 
is negligible and the depinning current density is given 
by the Eq. |S1 Let us also emphasize that the critical 
current of vortex depinning can not exceed the depairing 
one. Thus, Eqs. ([U E]) are valid only for jd < jc, i.e., 
for rather small size of the defect. In the clean limit, 
e.g., this restriction on the size reads: a < ^/q/kp ^ 
^. For larger defect radii the depinning current density 
should saturate at the jc value. The resulting schematic 
dependencies of critical current densities vs the defect 
radius a are shown in Fig. 2] for clean (a) and dirty (b) 
limits. One can see that the unusual bound vortex state 
near the defect can be observed only for rather clean 
samples with Aqt/H ^ 1. In the dirty case for Aqt/H < 
1 both critical current densities Jl and jd coincide. The 
depinning current estimate in this limit can be also found 
following the approach developed in Refs. [1, Q for the 
pinning on impurities with small scattering cross section. 

To sum up, we have developed a microscopic descrip- 
tion of the two-stage scenario of vortex depinning from 
small size columnar defect and predict the formation of 
stable vortex - antivortex configurations bound to the 
defect at intermediate transport current densities. The 
lower critical current density corresponds to the forma- 
tion of vortex -antivortex state and should determine 
rf nonlinear response of vortex system at rather small 
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transport current amplitudes. The upper critical current 
density jd corresponds to the destruction of the bound 
vortex configuration around the defect and should deter- 
mine the dcpinning transition at dc currents. 
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